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Summary 

New  mathematical  methods  have  been  developed  and 
prototyped  for  the  utilization  of  data  from  dynamic 
wind  tunnel  tests.  These  methods  provide  a  much 
improved  analytical  toolset  for  flight  dynamicists  to 
better  understand  aircraft  motion  and  stability  under 
extreme  maneuverability  conditions.  This  approach  is 
based  principally  on  the  technical  work  in  this  area  by 
the  late  Juri  Kalviste.  The  dynamic  stability  problem 
is  characterized  thoroughly,  from  data  generation 
through  practical  application. 

The  Phase  I  research  in  this  SBIR  concentrated  on 
identifying  and  proofing  a  set  of  analytical  tools  that 
provide  visibility  into  the  dynamic  stability  regime. 
Although  several  technical  applications  exist  for  the 
resultant  dynamic  data,  the  end  focus  of  this  research 
project  is  to  create  a  credible  aerodynamic  database 
usable  in  a  high  fidelity  6DOF  simulation.  The 
accomplishment  of  that  goal  by  the  end  of  Phase  11 
will  have  considerably  furthered  the  technical 
knowledge  base  in  this  field. 

The  technical  approach  for  the  project  can  be  divided 
into  three  main  parts, 

(1)  Dynamic  Data  Generation  Methods 

(2)  Data  Reduction  Methods 

(3)  Aerodynamic  Modeling  Methods 

The  core  of  this  SBIR  is  the  Data  Reduction  Methods. 
Necessary  attention  is  given  to  Dynamic  Data 
Generation  (the  source  of  the  data  to  be  reduced)  and 
Aerodynamic  Modeling  (the  application  of  the 
reduced  data)  in  order  to  provide  a  global  treatment  of 
the  end-to-end  problem. 

A  computer-aided-engineering  (CAE)  software 
program  has  been  prototyped  in  Phase  I  making  use  of 
the  reduction  algorithms  developed.  Once  completed 
and  commercialized,  the  product  name  will  be 
AeroMaster,  A  new  model-based  library  architecture 
has  been  applied  to  the  program.  Each  of  the  three 
areas  can  be  thought  of  as  an  independent  software 
model.  The  models  share  data  through  well-defined 
data  interfaces,  or  model  specifications.  This 
architecture  provides  a  very  powerful  technical 
advantage.  It  allows  for  different  models,  using 
different  methods,  to  be  directly  integrated  without 
affecting  the  overall  operation  of  the  program. 


For  example,  a  new  data  reduction  scheme  can  be 
programmed  as  a  standalone  model,  then  integrated 
into  AeroMaster  and  immediately  evaluated.  From  a 
User’s  standpoint,  the  tool  provides  the  best  of 
adaptability  and  expandability. 

Dynamic  Data  Generation  Methods 

The  dynamic  stability  parameters  consist  of  three 
rotation  rate  effects  (body  rates  P-Q-R),  and  two 
aerodynamic  angle  rate  effects  (a',  P'). 

Dynamic  stability  data  can  be  generated  either 
empirically  or  analytically.  The  empirical  methods 
include  tunnel  testing,  flight  test,  and  the  use  of 
analogous  data  (from  similar  configuration  vehicles). 
The  analytical  methods  consist  of  various 
computational  flow  models  that  consider  rotational 
and  accelerated  motion. 

The  Data  Reduction  model  is  designed  to  accept  time- 
history  state  data  from  any  data  generation  source. 

Data  Reduction  Methods 

Generalized  reduction  methods  have  been  developed 
to  provide  improved  signal  extraction  from  dynamic 
test  data.  The  models  use  curvefit  algorithms  tailored 
to  the  types  of  motion  used  in  rotary  balance  and 
forced  oscillation  tests.  The  two  methods  presented 
employ  a  linear  regression  fit  and  a  Fourier-series  fit. 
Both  models  incorporate  smoothing  and  quality-of-fit 
features.  These  methods  allow  the  practicing 
dynamicist  to  better  evaluate  test  data  and  correctly 
choose  the  data  format  for  an  aerodynamic  database. 

Aerodynamic  Modeling  Methods 

A  generalized  aerodynamic  database  and  coefficient 
buildup  method  is  presented,  suitable  for  both  linear 
and  non-linear  aerodynamic  data.  For  dynamic 
stability  parameters,  the  model  incorporates  a 
generalized  vector  interpolation  scheme  known  as  the 
Kalviste  Method.  The  model  is  database-driven  and 
designed  for  rapid  reconfiguration. 
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Problem  Statement 


The  correct  modeling  of  vehicle  aerodynamics  is 
essential  for  three  reasons: 

(1)  Technical  understanding  of  vehicle  flight 
behavior  to  perpetuate  the  design  of  more 
efficient,  safe,  and  maneuverable  aircraft. 

(2)  Training  of  aircrews  to  maximize  the 
flight  performance  of  a  given  vehicle  design. 

This  is  especially  important  in  simulator 
training  where  the  core  vehicle  simulation 
should  very  closely  match  the  actual  vehicle 
behavior  throughout  the  design  operating 
envelope, 

(3)  Design  of  flight  control  systems  that 
maintain  adequate  control  throughout  the 
operating  envelope  without  limiting  the  full 
performance  capability  of  the  vehicle. 

It  is  well  understood  that  actual  vehicle  aerodynamics 
are  influenced  by  the  dynamic  effects  of  rotational 
motion  and  translational  acceleration.  In  terms  of 
aerodynamic  stability  parameters,  these  dynamics  are 
modeled  as  effects  of  P-Q-R  and  a',  P'.  The 
quantitative  evaluation  of  the  dynamic  motion  terms 
has  historically  been  based  on  special  wind  tunnel 
tests.  These  tests  include  rotary  balance,  traditionally 
used  for  spin  analysis,  forced  oscillation,  plunging 
oscillation,  and  combined  motion  tests. 

The  dynamic  data  from  these  tests,  when  used  to 
predict  the  motion  of  an  aircraft  maneuvering  at  high 
angles-of-attack,  has  failed  to  correlate  well  with 
recorded  flight  test  behavior.  As  a  result,  the  6DOF 
simulations  used  by  major  aerospace  firms  for  flight 
analysis  and  aircrew  training  do  not  largely 
incorporate  the  tunnel  data.  The  special  apparatus  and 
test  procedures  for  dynamic  testing  is  costly,  so  one 
might  question  the  continuation  of  such  testing,  or  at 
least  the  continued  use  of  unsuccessful  analysis 
methods. 

The  dilemma  is  that  flight  test  provides  a  very  scant 
source  of  stability  data,  whereas  tunnel  test  provides  a 
vast  amount  of  stability  data.  If  simulations  are  going 
to  accurately  model  dynamic  effects,  the  data  is  likely 
to  come  from  tunnel  test.  Additionally,  from  an 
engineering  intuition  standpoint,  dynamic  tunnel  tests 
appear  to  be  a  valid  approach. 


From  a  macro  viewpoint,  the  deficiencies  could  be 
with  any  of  the  three  areas  previously  mentioned  - 
tunnel  tests,  test  data  reduction,  or  aerodynamic 
modeling.  On  the  basis  of  data  from  rotary  balance 
tests,  Juri  Kalviste  showed  that  the  dynamic  stability 
parameters  from  forced  oscillation  tests  could  not  be 
treated  as  linear  derivatives  at  flight  conditions 
involving  high  angles  of  attack.  He  also  showed  that 
that  aerodynamic  modeling  of  tunnel  data  produces 
vastly  different  (and  incorrect)  coefficients  when 
summed  using  linear  interpolation  vs.  vector 
interpolation. 

The  solution  to  many  difficult  technical  issues  lies  in 
having  a  fundamental  understanding  of  the  physics  of 
the  problem.  Many  engineering  and  scientific 
disciplines  have  successfully  used  scientific 
visualization  as  a  method  of  providing  that 
fundamental  insight.  The  idea  is  that  given  a  sufficient 
visual  model  of  the  technical  data,  the  mind  can 
heuristically  fathom  and  solve  the  problem. 

The  AeroMaster  software  program  uses  that 
philosophy  as  the  basis  for  its  design.  The  research 
performed  under  this  SBIR  is  developing  improved 
math  methods  for  reduction  and  aerodynamic 
modeling,  but  the  software  is  architected  to  accept 
alternative  methods  as  they  become  available  from 
government  or  industry  engineering  organizations.  It 
is  believed  that  this  tool  will  provide  the  necessary 
visibility  to  understand  and  solve  the  dynamic 
modeling  issue. 
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1.0  Wind  Tunnel  Test  Methods 

The  aerodynamic  characteristics  of  a  vehicle  can  be 
described  by  six  nondimensional  aerodynamic 
coefficients.  These  consist  of  three  force  coefficients 
and  three  moment  coefficients^  all  referenced  to  the 
principal  body  axes  X-Y-Z. 

Wind  tunnel  tests  are  performed  to  measure  the 
aerodynamics  of  a  given  air  vehicle  design  operating  at 
given  flight  conditions.  The  vehicle  models  are 
supported  by  a  sting  or  strut  apparatus,  and  force- 
moment  measurements  are  made  by  mass  balance  or 
strain  gauge  sensors.  The  dimensional  measurements 
are  converted  to  the  more  useful  and  compact 
nondimensional  coefficient  form  for  development  of 
an  aerodynamic  database. 

Tunnel  tests  can  be  grouped  into  two  categories,  static 
tests  and  dynamic  tests.  For  quasi-steady-state  motion, 
the  aerodynamic  coefficients  are  functions  of  three 
effects:  flight  condition,  aircraft  configuration,  and 
dynamic  motion  parameters. 

The  flight  condition  is  a  function  of  angle-of-attack 
(a),  angle-of-sideslip  (p),  Mach,  and  altitude. 

The  aircraft  configuration  is  a  function  of  a  given 
combination  of  control  surface  deflections,  wing 
sweep  position,  landing  gear/flap  positions,  etc. 

The  dynamic  motion  parameters  represent  body 
rotation  rates  and  aerodynamic  angular  rates.  The 
three  rotational  rates  about  the  principal  body  axes  are 
Roll  Rate  P,  Pitch  Rate  Q,  and  Yaw  Rate  R.  The  two 
aerodynamic  angle  rates  are  Angle-of-Attack  Rate  a' 
and  Angle-of-Sideslip  Rate  P'. 

Static  Tests 

Static  tests  are  a  function  of  flight  condition  and 
aircraft  configuration.  As  the  name  implies,  the 
aerodynamics  are  measured  with  the  model  in  a 
stationary  condition.  The  test  vehicle  is  configured 
with  a  given  geometry  and  positioned  at  a 
predetermined  angle-of-attack  and  angle-of-sideslip. 
Tunnel  flow  conditions  determine  the  test  Mach  and 
altitude. 


Dynamic  Tests 

Dynamic  tests  are  performed  to  measure  vehicle 
characteristics  while  undergoing  motion,  typically  for 
determining  spin  or  maneuvering  behavior.  The  two 
common  dynamic  tests  currently  in  practice  involve 
either  steady-state  or  oscillatory  motion.  Rotary 
Balance  tests  are  steady-state  in  nature,  while  Forced 
Oscillation  and  Plunging  Oscillation  tests  are 
oscillatory  in  nature. 

In  a  Rotary  Balance  test,  the  model  is  rotated  at  a 
constant  rate  about  an  axis  parallel  to  the  freestream 
velocity  vector  of  the  wind  tunnel.  At  the  reference 
body  origin  (RBO),  the  model  aerodynamic  angles  (a, 
P)  are  constant  with  respect  to  the  airstream 
throughout  a  rotational  cycle.  This  steady-state 
rotation  is  the  wind  axis  roll  rate.  Normally  this  data 
is  used  for  spin  analysis  since  the  rotary  motion 
closely  approximates  actual  spin  motion.  The  data  can 
also  be  used  for  maneuvering  analysis.  Care  must  be 
taken  to  ensure  that  it  is  used  in  the  same  manner  in 
which  it  was  measured.  A  vehicle’s  total  rotation 
vector  consists  of  components  that  may  include  both 
steady-state  rotation  and  oscillatory  terms  from  body 
PQR  rates.  Rotary  Balance  data  applies  only  to  the 
velocity  vector  component  of  the  total  rotation 
vector. 

In  a  Forced  Oscillation  test,  the  model  is  positioned  at  a 
predetermined  angle-of-attack,  then  forced  to  oscillate 
about  a  principal  vehicle  axis.  The  frequency  and 
amplitude  of  the  oscillation  determines  the  maximum 
rotation  rate  (P-Q-R).  This  rotation  about  a  body  axis 
also  causes  a  change  in  angle-of-attack  (a)  and  angle-of- 
sideslip  (P),  which  results  in  aerodynamic  angle  rate 
effects  (a',  P').  These  aerodynamic  angle  rate  terms 
have  both  rotational  and  translational  parts.  Only  the 
rotational  part  is  active  in  a  forced  oscillation  test. 

In  a  Plunging  Oscillation  test,  the  model  is  positioned 
at  a  predetermined  angle-of-attack,  then  forced  to 
oscillate  along  a  wind  axis.  This  measures  the 
translational  part  of  the  aero  angle  rate  terms. 

The  equations  for  static  and  dynamic  parameters  are 
provided  in  the  Reference  Equations  supplement  to 
Section  1.0. 
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Preliminary 


Reference  Equations 

1.0  Wind  Tunnel  Test  Methods 


Preliminary 


Rotary  Balance  Test 


Vj  =  Vdocity  Vector 

P  =  Rdl Rate  abmt  Body  Y-axis 
Q  =  Pitch  Rate  almt  Body  Y-axis 
R  =  Yaw  Rate  aboat  Body  Y-ax  is 

CO  =  Rotation  Rate  about  Vdocity  Vector 
CC  =  An^cf-Attack 
P  =  An^cfSideslip 


(a>)  =f(a,fi,P,qjt) 

0)  =  Poosaoosfi  +  Clsm/3  +  RsinaoosP 
0)  ^  (P^+  Q2+ 

(P,Q,R)  =f  (a),  a,  p) 

P  =  0)  axa  oosP 
Q  =  0)  sinP 
R  =  (O  sina  oosP 
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Preliminary 


Reference  Eaucttions 

1.0  Wind  Tunnel  Test  Methods 


Preliminary 


Forced  Oscillation  Test  about  Body  X-axis 


X 


B 


Yn 

Yb 


=  Body  X-axis 
=  Body  Y-axis 
—  Body  X-axis 


(X  =  A  n^cf-A  ttack 
P  —  An^cfSideslip 


B 


6)  =  OsdlatimFvaqtioicy 
(j)  =  EderRoHAn^ 

P  —  RoU  Rate  about  Body  Y-axis 


(Pmax)  =  /H 

P MAX  =  ®  *I^MAX 


(at)  -  f  (Pmax>  ^max) 
^  =  Pmax  Y  ^max 


'yP)  —  f  ft)  ^MAx) 

<!>  =  (I>max^^(o>  t) 

<t>'  =  P  ’-^MAXOi  COs((Ot) 

^«  =  F  =-(l)MAx(0^^in(G)t) 


(a,p)  =f(ao,^) 

tana  =  tanUQCOS^ 
sinfi  =  smaQsiru^ 


(a-.fi-)  =f(a,fi,P) 

a‘  =  ■  P  oosa  tanfi 


P'  =  Psina^ 
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2.0  Data  Measurement  and  Reduction 

Measurements  in  a  wind  tunnel  are  performed  with 
sensors  mounted  in  the  vehicle  support  apparatus,  the 
sensing  device  generally  being  a  strain  gauge.  The 
sensors  are  mounted  in  pairs  with  known  moment 
arm  geometry,  thereby  enabling  measurement  of  the 
forces  and  moments  experienced  by  the  test  vehicle. 
The  data  are  measured  at  some  predetermined  sample 
rate  over  some  sample  interval.  Signal  processing  is 
performed  on  the  raw  data  for  analog-to-digital  (ADC) 
and  engineering  unit  (EUC)  conversions. 

As  with  all  measurements,  inaccuracies  arise 
throughout  the  process  in  the  form  of  noise,  bias, 
histeresis,  etc.  The  reconstruction  of  signals  into 
meaningful  data  requires  an  understanding  of  the 
physics  of  the  waveform  and  the  source  and 
magnitude  of  measurement  errors.  The  selection  of 
appropriate  analysis  techniques  can  provide  natural 
correction  and  smoothing  of  raw  data  to  better  extract 
the  inherent  signal. 

For  the  generation  of  an  aerodynamic  database,  the 
processed  parameters  are  stored  in  tables  that  are 
organized  in  a  manner  that  is  both  readable  and 
conducive  to  rapid  search  and  interpolation.  The 
conventional  method  of  database  organization  stores 
non-dimensional  force-moment  coefficients  in  blocks 
of  baseline  coefficients,  supplemented  by  blocks  of 
incremental  coefficients.  The  baseline  coefficients 
represent  static  data  for  a  given  flight  condition.  The 
incremental  coefficients  represent  static  data  for  a 
given  flight  condition  and  aircraft  configuration^  or 
dynamic  data  for  a  given  flight  condition^  aircraft 
configuration  and  dynamic  motion. 

For  static  and  dynamic  wind  tunnel  tests,  two 
techniques  are  presented  to  provide  smoothing  and 
signal  extraction  on  measured  data.  For  non-cyclical 
tests  (^tatic  and  rotary  balance),  a  linear  regression 
algorithm  is  applied  due  to  the  steady-state  nature  of 
the  data.  For  cyclical  tests  {forced  oscillation  and 
plunging  oscillation),  a  Fourier  series  algorithm  is 
applied  due  to  the  oscillatory  nature  of  the  data. 


Linear  Regression  Algorithm 

Regression  methods  are  used  to  curvefit  sampled  data 
that  is  curvilinear  ascending  or  descending.  The  four 
classic  regression  curvefits  are. 


Straight  Line: 
Exponential  Curve: 
Logarithmic  Curve: 
Power  Curve: 


y  ^  a  +  bx 
y  =  ae^ ,  (a>0) 
y  ^  a  bln  X 
y  =  ax^,  (a>0) 


All  provide  convenient  equations  to  compute  values 
for  the  dependent  variable,  y,  given  any  value  of  the 
independent  variable,  x. 


The  linear  regression  fit  (Straight  Line)  is  applicable  to 
non-  cyclical  wind  tunnel  tests.  When  implemented  in 
an  automated  computer  program,  the  method  solves 
for  the  coefficients  (a,  b),  and  provides  a  quality  of  fit 
coefficient  (r^). 

The  values  of  the  coefficients  are  indicative  of  the 
cleanliness  of  the  sampled  data.  Ideally,  the  coefficient 
b  should  be  near  zero,  indicating  the  steady-state 
nature  of  the  data.  The  coefficient  of  determination  r^ 
should  be  close  to  1.00,  indicating  a  very  good  fit. 
The  coefficient  a  is  the  desired  steady  state  value  of  the 
measured  parameter. 

The  reduction  process  flow  is  provided  below. 

(1)  Enter  n  samples  of  ( x ,  y )  data  pairs 

(2)  Compute  the  coefficient  b 

(3)  Compute  the  coefficient  a 

(4)  Compute  the  coefficient  r^ 

(5)  For  Rotary  Balance  data,  compute  the 
incremental  force-moment  value  by 
subtracting  the  static  value  at  the  same 
flight  condition. 

(6)  Compute  the  normalized  force-moment 
coefficient  value. 

(7)  Store  the  coefficient  in  the  aero  database  table. 

The  linear  regression  equations  are  provided  in  the 
Reference  Equations  supplement  to  Section  2. 


Fourier  Series  Algorithm 

Fourier  methods  are  used  to  curvefit  sampled  data  that 
is  cyclical  in  nature.  For  forced  oscillation  tests  where 
the  attributes  of  the  oscillatory  motion  are  known,  a 
technique  has  been  developed  to  fit  the  sampled  data 
in  a  manner  that  provides  good  signal  extraction  and 
natural  smoothing. 
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The  objective  of  a  forced  oscillation  test  is  to  measure 
the  effects  of  rotational  or  translational  motion  on  the 
aerodynamic  forces  and  moments  experienced  by  the 
aircraft.  Two  technical  issues  require  attention.  The 
first  is  the  extraction  of  the  desired  parameter  from 
the  sampled  data.  The  second  is  the  analysis  of  the 
parameter  to  determine  whether  the  effect  of  motion 
is  near-linear  or  clearly  non-linear.  This  latter 
determination  will  decide  how  the  data  will  be  stored 
and  retrieved  in  an  aerodynamic  database. 

A  tailored  Fourier  analysis  method  addresses  the  signal 
extraction  and  smoothing  issue. 

Given  an  array  of  raw  time-sampled  data,  a  sample  set 
N  is  chosen,  where  N  is  the  number  of  data  points  per 
cycle.  Ideally,  the  value  of  N  should  be  as  no  greater 
than  2x  the  highest  expected  frequency  of  the  native 
signal  (per  Nyquist  theory).  Given  some  knowledge 
of  the  expected  noise  amplitude  A,  and  desired  noise 
error  E,  the  number  of  cycles  k  may  be  computed  to 
achieve  the  desired  signal-to-noise  (S/N)  value.  A 
Fourier  series  approximation  is  applied  to  the  sample 
set  and  the  Fourier  coefficients  extracted. 

In  equation  form, 

L  Oscillation  Period  SEC 

k  Oscillation  Cycles  Cycles 

y  Signal  Amplitude 

A  Noise  Amplitude 

E  Noise  Error 

N  No.  of  Pts  per  Cycle 

n  Sample  Index 

t  Time  Point  SEC 

Fourier  Coefficients 

Number  of  sample  cycles  determination, 

k  =A/{N^E) 

Fourier  series  approximation. 

Trigonometric  form, 

F(t)  =  sin( 2mt/L)  -i-  cos( 2mt/L) 

Exponential  form, 

F(t)  = 


Once  the  Fourier  coefficients  are  computed  from  the 
raw  sample  data,  the  fit  is  complete.  The  frequency 
content  and  noise  elimination  were  accommodated  by 
the  selection  of  N  and  k. 

Much  can  be  extracted  from  the  final  Fourier  fit.  The 
value  of  the  measured  signal  at  any  time  t  can  be 
directly  computed.  For  a  force  or  moment  parameter, 
this  represents  the  total  value.  Since  the  Fourier 
term  is  a  constant  that  represents  the  value  with  no 
motion  effects,  this  is  equivalent  to  the  static  part  of 
the  parameter  value.  The  remaining  part  represents  the 
incremental  effect  due  to  motion,  which  is  appropriate 
for  storage  in  the  aerodynamic  database.  Since  the 
Fourier  fit  for  any  given  parameter  has  time  as  the 
independent  variable,  the  incremental  term  must  be 
correlated  with  the  appropriate  dynamic  parameter  to 
represent  an  aerodynamic  term.  For  instance,  in  a 
forced  oscillation  test  about  the  X-axis,  the  rolling 
moment  Q  is  measured  as  a  function  of  time,  as  is  the 
roll  rate  P.  The  desired  aerodynamic  term  is  ‘rolling 
moment  as  a  function  of  roll  rate’,  or  Q  (P)y  requiring 
a  correlation  of  Q  (t)  and  P(t). 

The  reduction  process  flow  is  provided  below. 

(1)  Choose  a  value  of  N  based  on  motion  rates. 

(2)  Compute  the  number  of  sample  cycles  k, 

(3)  Enter  N samples  o{{y  ,t)  data  pairs  per  cycle 
at  a  sample  interval  of  L/ N  over  the  raw  data, 

(4)  Compute  the  Fourier  coefficients  and  b^ 
over  k  cycles. 

(5)  Subtract  the  term  (static  part). 

(6)  For  motion  parameters,  determine  rate  and 
acceleration  values. 

(7)  Correlate  the  force-moment  term  with 
the  appropriate  dynamic  motion  term. 

(8)  Generate  values  for  the  aerodynamic  database. 
For  near-linear  data  generate  a  derivative. 

For  non-linear  data,  generate  an  ascending  data 
table  of  the  aero  term  vs.  the  motion  term. 

(9)  Compute  the  normalized  force-moment 
coefficient  value. 

(10)  Store  the  coefficient  in  the  aero  database  table. 

The  Fourier  equations  are  provided  in  the  Reference 
Equations  supplement  to  Section  2.0. 
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Preliminary 


Reference  Equations 

2.0  Data  Measurement  and  Reduction 


Preliminary 


Linear  Regression  Fit 


X  =  Inkperdent  VariaUe 
y  =  Dependent  VamUe 
Ji  ”  iV 


a  =  Y-okis  Interoept 
h  =  Slope 

=  Co^lkient(fDetemimtk^ 


r2  = 

„  IX:  Sy: 

exiji  ■  ‘ 

2 

X  2  .  (^Xi)  ^ 

^  n 

ry2  .  (^yi)^ 
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Preliminary 


rt  r  T’  •  Preliminary 

Reference  Equations 

2.0  Data  Measurement  and  Reduction 


-  Fourier  Series  Fit  - 

i< -  ^  - >1 


TrigjmretricForm  Exponent  Form 


a„  <=  1/N  Z  y(t)  ox(2nrt/ L) 

y(t)  -Ed.- 

=  1/N Z  y(t)  sm(2nrt/ L) 

y(t)  =  Zd„  ^(-2n  ^j^n) 

f(t)  =Zd„  -(-271  -j-nf  - 

sm(/>„  =-b„/(a„^+ 

=  FcMnerCo^tdent 

y(t)  =  Z  a,^cos  (27tnt/ L)  +b^sin(2nnt/ L) 

n  =  Sanplelnchc,  0  toN-1 

]  =  Irm^mryNimher 

y(t)  =  Z  c„ocs  (2rrnt/L  +  <f)J 

y(t)  =  ValuecfDV  attirrEt 

Frequency  Phase 

y(t)  =  FmtDerimtkeatUmet 

^n>  K>  K  =  Fourier Co^Jkiencs 

y^(t)  =  Second  Derkatke  at  Urne  t 

y(t)  =  Value  (fDV  at  tmet 
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3,0  Aerodynamic  Data  Representation 

The  aerodynamic  characteristics  of  a  vehicle  can  be 
described  by  six  nondimensional  aerodynamic 
coefficients: 

Three  force  coefficients  Cx,  Cy,  Cz 
Three  moment  coefficients  Q,  Cj^,  C„ 

Each  total  coefficient  can  be  partitioned  into  two 
parts,  a  static  coefficient  and  a  dynamic  coefficient. 

The  static  coefficients  are  a  function  of  flight  condition 
and  aircraft  configuration^  and  represent  non-rotational 
parameters. 

The  dynamic  coefficients  are  a  function  of  flight 
condition^  aircraft  configuration^  and  dynamic  motion 
parameters. 

For  any  of  the  six  coefficients, 

CtOTAL  =  ^STATIC  ^DYNAMIC 

where 

^STATIC  “  /  ( oc,  P,  M,  h,  5 ) 

^dynamic  /^  ( Pj  M,  h,  5,  P,  Q,  R,  06 ,  p  ) 

The  generalized  form  is. 

Coefficient  Variable  ==  /  (Independent  Variables) 
or 

CV  =/(IV’s) 

In  this  form  the  coefficients  are  represented  as  a  single 
non-linear  function  of  all  of  the  variables.  Although 
this  would  result  in  a  very  accurate  aerodynamic 
model,  it  presents  an  impractical  solution  due  to  the 
enormity  of  the  resulting  database.  A  more 
convenient  method  defines  each  total  coefficient  as  a 
sum  of  component  terms.  This  is  known  as 
coefficient  buildup. 

Two  approximations  are  made  to  simplify  the 
representation  of  aerodynamic  data.  The  first 
approximation  is  to  separate  the  total  CV  into 
incremental  coefficient  terms  that  sum  to  the  total, 
where 

CVtotal  ==  ACVi  +  ACVz  +  ...  ACV^ 

The  second  approximations  is  to  treat  the  non-linear 
incremental  terms  as  products  of  a  linear  derivative 
and  an  independent  state  variable,  where 


ACVi  =  /  (Dependent  Variable  x  Independent 
Variable) 
or 

ACVi=/(DV  xIV) 

Care  must  be  taken  with  both  approximations.  The 
first  assumes  that  the  Independent  Variable  (FV)  of 
each  incremental  coefficient  term  is  unaffected  by  the 
state  change  of  another  IV.  The  second 
approximation  assumes  that  the  vehicle’s 
aerodynamics  for  a  given  ACV  are  a  linear  function  of 
its  corresponding  IV.  The  decision  as  to  which 
simplifying  approximations  are  tolerable  is  based  on 
the  magnitude  and  nonlinearity  of  each  DV.  This  can 
often  be  determined  with  good  computational 
methods,  but  most  commonly  by  wind  tunnel  testing. 

Standard  Aerodynamic  Model 

A  generalized  method  has  been  developed  that  allows 
coefficient  buildup  to  be  done  using  database  entry  in 
lieu  of  software  coded  equations.  Maximum  flexibility 
is  provided  in  the  coefficient  buildup  process,  and  the 
result  is  an  intuitively  organized  database. 

Each  total  coefficient  is  computed  using  a  “sum  of 
products”  method,  or  mathematically, 


Each  incremental  coefficient  term  may  represent  a 
constant,  a  linear  derivative  term,  or  a  non-linear 
term.  The  symbology  for  each  option  is, 

ACVi  =  K  Constant  Term 

ACVi  =  DViv  X  rV  Linear  Derivative  Term 
ACVi  ^  DV  (rV’s)  Nonlinear  Term 

The  buildup  of  drag  coefficient  is  provided  as  a  simple 
example,  where  the  total  coefficient  is  a  function  of 
three  incremental  terms,  each  represented  as  one  of  the 
three  ACV  options, 

CD  =  CDq  Parasite  Drag  Constant  Term 

+  aCDm  X  M  Mach  Drag  Derivative  Term 
-1- aCd  (a)  Induced  Drag  Nonlinear  Term 
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The  nonlinear  terms  can  of  course  be  functions  of 
more  than  one  independent  variable.  The  common 
method  of  storing  the  DV  data  is  in  large  tabular 
databases.  The  instantaneous  DV  value  is  obtained  by 
performing  linear  interpolation  based  on  the  state 
value  of  the  IV.  For  multiple  IV's,  a  multidimensional 
linear  interpolation  is  performed.  On  modern 
aircraft,  four-  or  five-dimensional  table  lookups  are 
common. 

Kalviste  Model  for  Dynamic  Data 

It  was  cited  in  Section  2.0  that  data  should  only  be 
used  in  accordance  with  how  it  was  measured. 
Applying  this  principle  to  measured  dynamic  data,  a 
method  has  been  developed  to  correctly  use  the  results 
of  rotary  balance,  forced  oscillation,  and  plunging 
oscillation  tests  in  a  single  aerodynamic  database. 

This  technique,  known  as  the  Kalviste  Method,  has 
evolved  since  the  early  1980’s.  The  reduced  method 
involves  resolving  the  total  rotation  vector  into  two 
components,  always  in  the  X-Y  plane.  The 
generalized  method  resolves  the  total  rotation  vector 
into  three  components,  and  considers  all  three  body 
axes.  One  component  is  the  steady-state  rotation 
along  the  velocity  vector  (cOmod)^  while  the  other 
components  are  the  corrected  body  rotation  rates 
(Pmodj  Qmod>  Rmod)*  ^  vector  addition  technique  is 
used  to  combine  coefficient  terms  from  different  W.T. 
test  methods,  and  is  key  to  the  overall  correctness  of 
the  aerodynamic  coefficients.  Where  the  velocity 
vector  and  rotation  vector  have  large  separation,  a 
difference  in  coefficient  buildup  techniques  can  result 
in  a  spin  vs.  non-spin  result. 

The  total  rotation  vector  is, 

|Q|  =  (P2  +  Q2  + 

The  generalized  three-dimensional  components  are, 

^MOD  =  P  ”  ®M0D  cosot  COSp 

Qmod  =  Q  -•  C0j^o£)  sinp 

P-MOD  —  R  -  0)mod  sina  cosp 

The  method  tests  all  combinations  of  rotation  rate 
magnitude  to  determine  which  three  components  to 
use.  One  or  more  of  the  four  terms  will  equal  zero. 


Once  the  elimination  procedure  is  complete,  the 
rotation  rate  (--)mod  parameters  are  treated  as  normal 
state  variables.  Entry  is  made  into  the  aerodynamic 
model  and  corresponding  coefficient  terms  are 
computed. 

To  illustrate  the  use  of  wind  tunnel  test  data,  a 
complete  coefficient  buildup  will  be  described  that 
uses  data  from  all  dynamic  tests.  This  method  is 
suitable  for  use  in  a  6DOF  simulation,  and  considers 
all  of  the  rotational  motion  effects.  The  buildup  uses 
the  approximation  that  the  total  coefficient  can  be 
partitioned  into  incremental  parts. 


-  Qbase 

Static  Test 

+  A  Ci  (co  j^od)  Rotary  Balance  Test 
+  A  Ci  (P  mod)  Forced  Oscillation  Test  (X-axis) 

+  A  Ci  (Q  mod)  Forced  Oscillation  Test  (V-axis) 

+  A  Ci  (Rmod)  Forced  Oscillation  Test  (Z-axis) 

+  A  Q  (a  ’  t) 

Plunging  Oscillation  Test  (Zj^axis) 

+  aQ  ((3‘t) 

Plunging  Oscillation  Test  (V ^axis) 

a  \ 

a  ‘  due  to  Translation 

P‘t 

P  ‘  due  to  Translation 

The  aero  rate  terms  (a  P  ')  have  both  rotational  and 
translational  parts.  By  the  physics  of  the  test  motion, 
the  forced  oscillation  data  include  the  rotational  part  of 
the  aero  terms  r  ,  r  ).  The  complete  coefficient 
buildup  requires  inclusion  of  the  remaining 
translational  part.  These  terms  can  be  measured  in  a 
translational  acceleration  tunnel  test,  or  plunging 
oscillation  test,  where  the  aircraft  attitude  is  fixed  and 
the  velocity  vector  oscillates. 

The  expanded  coefficient  equations  are  provided  in  the 
Reference  Equations  supplement  to  Section  3.0. 
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Preliminary 


Reference  Equations 

3.0  Aerodynamic  Data  Representation 


Preliminary 


Aerodynamic  Coefficient  Model 


P  P 

MOD  ^  ^ 


R  Y  ^MOD 


(Total  Rotation  Vector) 


Q  =  Total  Rotation  Vector 
CO^  —  Rotation  Rate  aboat  Vdocity  Vector 
Pmod  Rdl  Rate  about  B(xfy X-axis  (Modftec^ 

Qmcd  “  Pitdo  Rate  about  Baiy  Y-axis  (Mallfie^ 
Rmod  ~  YawRate about  Body  2-axis  (Modfied) 


V^(  vdocity  Vector) 


Coefficient  Buildup  with  Dynamic  Test  Terms 
C  =  C 

BASE 

+  /iQ  (  CO^cD  Vr)  Rotary  Balance  Test 

+  Aq(P^^b/2Vr)  ' 

+  aQ  ( Qmod  c/2Vy)  •  TcnedQdllationTest 

+  (Rmod  t/2Vr)  , 


+  Aq(  aj‘  C/2V-J') 

+  Aq(  Pt'  t/2Vr) 


Plun^ng  Oscillation  Test 


b  =  R^erenoe  Lengh  along  Span 

c  =  R^erenceLengh  along  Oxmi 

Vj  =  True  Airspeed 
Op  *  y  Pt'  ^  Translatioal  part  cf  a\  p\ 

Note  1 :  Forosd  Oscillation  Test  terms  indude 
rotatioal  effects  ap\  pp\ 

Nctel:  Any (o^fcient terms fro^ 

substituted  mth  oonpuiational  terms. 


Kaluste  Rotation  Rates 
^MOD  ^  ^  “  ^MOD 

Qmcd  =  Q  ‘  ^moo  ^^'^P 

^MCD  =  ^  ■  (Omod^^'^^P 

Kahiste  Vector  Elimination  Procedure 

[Stepl]  =  P  /oosacosp 

Test:  sigi(q^^)  =sigi(Q)  and  sign  (Rj^^j^)  ^sigifR) 

Test:  Q^qj)  <Qand  Rj^^qd  <R 

Pass :  Use  Pj^q^  =  0 ,  ^mcd>  ^mod 

[Step  2]  =  Q/sinP 

Test:  sign  (P^^)  ^  sign  (P)  and  sign  (R^^)  ^sign(R) 

Test:  Pj^Q[)  <P and  Rj^qjj  <R 

Pass :  Use  =  0 ,  Pi^ao^  ^mqd 

[Steps]  (x>^^  =  R  /sinaoosfd 

Test:  sigz(P^^)  - sigi(P) andsigt(q^^j^)  =sigi(Q) 

Test:  Pj^^(jQ  <P and  <Q 

Pass :  Use  R^(^  ==  0  y  Pmo)^  Qm(X)>  ^mod 

[Step  4]  Allfaib  orallzerv: 

Rmod  =  Qmcd  =  ^mod  ^  ^mod  =  ^ 
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4.0  Sample  Reduction  Plots 

The  prototype  AeroMaster  software  program  was  used  to  generate  various  time-history  files.  The  data  are  meant  to  be 
representative  of  that  measured  in  static  and  dynamic  wind  tunnel  tests.  The  plots  overlay  both  clean  and  noisy  signal 
data  in  order  to  test  the  linear  regression  and  Fourier  reduction  methods. 

The  time-history  plots  are  provided  in  the  Reference  Plots  supplement  to  Section  4.0. 


5.0  AeroMaster  Software  Program 

As  of  the  completion  of  Phase  I,  the  AeroMaster  software  is  in  an  alpha  pre-release  state.  Informal  testing  has  been 
conducted  on  the  models,  with  the  emphasis  on  the  reduction  model.  The  completion  of  the  design  will  occur  in 
Phase  n,  as  well  as  the  formal  testing.  Ongoing  verification  will  occur  as  actual  tunnel  test  data  is  imported,  and  an 
aerodynamic  database  generated  and  simulated. 

The  User  is  encouraged  to  exercise  the  program  and  provide  written  feedback  regarding  features,  friendliness, 
anomalies,  errors,  and  overall  usability. 
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Preliminary 


Reference  Plot  Data 

4.0  Sample  Reduction  Plots 


Preliminary 


Linear  Regression  Fit  -  Static  Test 


A  linear  regression  algorithm  is  used  as  the  fit  method  for  non- oscillatory  data.  The  plot 
shown  below  is  an  example  of  static  data  from  the  AeroMaster  software  program.  The 
aerodynamic  data  were  generated  with  and  without  s^nal  noise.  The  curves  are  labelled 
‘StaticGlean’  and  ‘RAW  -  StaticNois)^^it’,  respectively.  The  regression  fit  file  is  labelled 
‘StaticNois}i^it’.  Ideally,  the  overplot  of  the  clean  signal  vs.  the  regression  fit  of  the  noisy 
signal  would  overlay  each  other.  In  this  example,  the  curvefit  shows  both  a  bias  and  a  slope, 
indicating  some  residual  error  due  to  the  extreme  noise  level  (S/N  of  1.0).  The  flow  and 
noise  parameters  are  listed  below. 


Angle-of-Attack 

10° 

Angle-of-Sideslip 

0° 

Airspeed 

0.5  Mach 

Static 

Motion 

Samples  per  Cycle 

128 

Signal-to-Noise 

1.0 

^AetoMastei  -  Standaid  Aero  Gen.  FI  BC^FullPanel 


£#8  £ditASPIFI(!S 


mm 

0 

-0.0124 

'0.00763284 

-0.00763284 

0  -0.8966 

-0.865047 

-0.865047 

0 

-1.855E7e-005 

0.0058 

-0.00758171 

-0.00763284 

5.11277C-005  -0.8611 

-0.865053 

-0.865047 

-6.31809e-t 

-3.71076C-005 

-0.0192 

-0.00753058 

-0.00763284 

0.000102255  -0.9014 

-0.86506 

-0.865047 

-1.26958e-I 

-5.56633C-005 

0.016 

-0.00747345 

-0.00763284 

0.000153383  -0.8449 

-0.865066 

-0.865047 

-1.901396-1 

-7.421 52C-005 

0.008 

-0.00742833 

-0.00763284 

0.000204511  -0.8854 

-0.865072 

-0.865047 

-2.5391  Be-l 

nmz  ; 

-9.27709C-005 

0.0011 

-0.0073772 

-0.00763284 

0.000255639  -0.8482 

-0.865079 

-0.865047 

-3.1 7097e-I 

mm  1 

-0.000111327 

-0.0444 

-0.00732607 

-0.00763284 

0.000306766  -0.8931 

-0.865085 

-0.865047 

-3.802786-1 

-0.000129679 

0.0275 

-0.00727494 

-0.00763284 

0.000357894  -0.8546 

-0.865091 

-0.865047 

-4.440556-1 

! 

-0.000148434 

-0.0392 

-0.00722381 

-0.00763284 

0.000409022  -0.8196 

'0.865098 

-0.865047 

-5.072366-1 

^  Close. 

. 

. . ... 

T_AeroCoefficIent_CY  -  StaticClean.FId 
IT  AeroCoeffIclenI  CY  RAW  -  StallcNolsvFlt. 


hi  J  in!i  !i  S!  >\\  III 


.  . . .  _ _ 


. 
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Preliminary 


Reference  Pht  Data 

4.0  Sample  Reduction  Plots 


Linear  Regression  Fit  -  Rotary  Balance  Test 


A  linear  regression  algorithm  is  used  as  the  fit  method  for  non-oscillatorydata.  The  plot 
shown  below  is  an  example  of  static  data  from  the  AeroMaster  software  program.  The 
aerodynamic  data  were  generated  with  and  without  signal  noise.  The  curves  are  labeled 
‘RBQean’  and  ‘RAW  -  RBNois)i^it’,  respectively.  The  regression  fit  file  is  labeled 
‘RBNois}^^it’.  Ideally,  the  overplot  of  the  clean  signal  vs.  the  regression  fit  of  the  noisy  signal 
would  overlay  each  other.  In  this  example,  the  cuivefit  shows  both  a  bias  and  a  slope, 
indicating  some  residual  error  due  to  noise.  The  flow  and  noise  parameters  are  listed  below. 
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Improved  VWnd  Tunnel  Data  Reduction  Procedure 


Preliminary 


Reference  Plot  Data 

4.0  Sample  Reduction  Plots 


Preliminary 


Fourier  Series  Fit 


A  Foxirier  series  algorithm  is  used  as  the  fit  method  for  oscillatory  data.  The  plot  shown 
below  is  an  example  of  forced  oscillation  data  from  the  AeroMaster  software  program.  The 
aerodynamic  data  were  generated  with  and  without  serial  noise.  The  curves  are  labeled 
TCXlean*  and  ‘RAW  -  FONois)^^it’,  respectively.  The  Fourier  fit  file  is  labeled  ‘FONois}Fit’. 
Ideally,  the  overplot  of  the  clean  signal  vs.  the  Fourier  fit  of  the  noisy  signal  would  overlay 
each  other.  That  is  the  case  in  this  example.  The  flow  and  noise  parameters  are  listed  below. 
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^pendix  A 


SBIR  Phase  I  Progress  Charts 


Monthly  progress  briefings  in  the  form  of  viewgraph 
presentations  were  provided  as  part  of  the  Phase  I 
effort.  These  charts  are  included  in  this  appendix  as 
part  of  the  overall  documentation  for  this  project. 


i^pendix  A 


Sight,  Soundy  and  Motion 


A-1 


SBIRN96-019 


Improved  Vflnd  Tunnel  Data  Reduction  Procedure 


Progress  Report  #  1 


Sighty  Sound,  and  Motion 


SBIRN96-019 

(proved  Wnd  Tunnel  Data  Reduction  Procedure 


O  O 

Q 


C 

o 

s 

cp  ON 

t-H 

5  si) 

Q  :z; 


<N 

0 

^  y 

^  It  § 

ON  ^ 

s^< 

O  <L)  _X 

Ph  -R 
d  ”ci  ’-i^ 

<u  <u 

S  I  O 

S 


I  K  ^ 

^|6|^ 

o  I  3  Q  ;< 

S  K  §  a  > 
a  s  R  12  § 
§. 

1  ?j  ^ 


325  East  HiUcrest  Dr.  #  100 
Thousand  Oaks,  CA  91360-5828 


SBIRN96-019 

Improved  Vflnd  Tunnel  Data  Reduction  Procedure 


C/5 

C/5 


□  on  □  □  □  □ 


Sight,  Sound,  and  Motion 


SBIRN96-019 

mproved  Wind  Tunnel  Data  Reduction  Procedure 


Sight,  Sound,  and  Motion 


Wind  Tunnel  Data  or  Water  Tunnel  Data.  The  data  can  be  static  or  dynamic. 
Dynamic  data  represent  the  increments  due  to  rotational  motion. 
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(2)  Design,  program,  and  test  a  math  model  to  corteedy  combine  contributive 
coefficient  terms  from  various  test  sources. 
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Plunging  Motion  Tests  (Translate  along  'Wind  axes) 

The  model  is  set  at  selected  a,  |3,  conditions  and  oscillated  at  a  selected  frequency-amplitude 
combination  along  a  selected  wind  axis. 
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Velocity  Vector  Component :  Q  mod  =  R  /  sin  a  Velocity  Vector  Component :  Q  mod  =  P  /  cos  a 

Body  X-axis  Component :  P  mod  =  P  -  Q  mod  (  cos  a  )  BodyZ-axis  Component :  R  mod  =  R  -  Q  mod  (  sin  a  ) 

Body  Z-axis  Component :  R  mod  =  0  Body  X-axis  Component :  P  mod  =  0 
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The  term  due  to  rate  of  change  of  Sideslip  ( p  trans  ),  provides  an  increment  for  the  translational  effects 
only.  The  rotational  effects  ( p  rot  )  are  included  in  the  rotary  term  (P  mod  ,  Rmod)  coefficients. 


a 

O 

a 

c 

o 

J 

0^  (U 

S 

SD  S 

za 

Pi  'S 

t-H  g 

^  a 

c/)  3 
H 

a 

TS 

S 

I 


a 

o 

o 

a 

o 

9J 

•c 

CJ 

CA 

<L> 

Q 

<u 

u 

Sfi 

<L> 

a 

a 


(U 


O 

00 


oj 


CL,  ^ 


•3  a 


42 

u 

p 

"p 

a 

PLh 


I 

c/5 


tu 

w 


03 


03 

+ 


<N 


03 

+ 

Co 


<N 

R. 

II 

•^A 

03 


U 


c 

I 

a 

"S 

0) 

u 

a> 

O 

u 

I 


03 


U  U 


•an. 

CS 

•  s 


a 

w  q; 

I 

G  a 


ca 


-  6 


1 

o 

I 

Ji 

M 

H 


a 

*13 

I 

a> 

Q 

II 


jj 

3 

«« 

•c 

"S 

s 

<D 

"B 


Jh  ^ 

^  M  s 


^  s:  iy 

■a^fS 


CJ 


^  a  § 

CD 


:r  X  ^  ^ 

Mh  U  1 — I 


Sight,  Sound,  and  Motion 


I 


0^ 

C 

Q 

S  S 
Z  Q 

HH  g 

^  I 

c/)  3 


C 

s 

I 


a 

o 

•^3 

o 

4 

n 

n 

Q 

4-> 

C/) 

q> 

H 


C/3 

G 

O 

•Z3 


V3 

a> 

H 


S 

H 


I 

cu 

Q 

§ 


CO 


•  ca 

t\ 

•  S 

a> 

ca  1=^ 

b  eC 


52 

1 

2 
CS 

Ph 


>  pL.  Ph 


Q 


52 

I 


(D 

CO 

os 


»H 

a 


52 

ii 

1 

2 

PP 

a> 

CA 

a 

(S 


§ 


52 

B 

1 

2 

v> 

a> 

H 


I 

CD 


CO 

<U 

H 


;:i 

H 


8 


CO, 


a 

u 

< 


a 

r\ 

CO. 

d 


a> 

CJ 

a 

’rt 

pq 

b 

B 

c§ 


8 


8 


8 

s 

Ph 

CN 

CSl 


a" 

<1 


pq 

c? 

CP 

rv 

ca 

r\ 

b 


c 

o 

•X3 


o 

CO 

o 

"T3 

(U 

fc5 

o 

Ph 


'  ca 


H 

•  d 

ca 


•  ca 

•b 

#v 

ca 


u 

< 


c\ 

ca 


a 

o 

‘S 


bD 

a 

*3) 

c 

3 


Sight,  Sound,  and  Motion 


SBIRN96-019 

Improved  Tunnel  Data  Reduction  Procedure 


a 

o 

*^3 


f/i 


(U 

H 


5h 

g 

M-h 

•I 

c/5 

r?5 

c/5 

<u 

<U 

P9 

b 

(2 

a 


jj 

CTS 

i 

■f . 

s  a 

C  <U 
CL,  JH 

I -I 

<U  ^ 
c/5  y 
CTS  Ch 

^  r> 
tw  .-y 

m-h 
^  (U 

•l  ° 

ll 

1 1 

<u  a 

ctS  6 
(U  O 


cu 


o 


a 

o 


■B 


aQ  =  f(P) 

[4]  Generate  a  database  of  the  Coefficient  (DV)  as  a  The  independent  state  variables  are  a,  R,  and  P.  For  the 

function  of  the  Flight  Orientation  and  Dynamic  between  the  +X-axis  and  the  velocrty  vector, 

Parameter  (IV).  (See  Example)  the  database  of  is  entered  with  a  and  P  mod.  Interpolation  is 

performed  to  determine  A  C  . 
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Translational  Aero  Rates.  The  rate  of  change  of  a  or  P  due  to  vehicle  translation,  due  to 
a  change  in  direction  of  flight. 
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Sum  of  all  contributing  rotation  terms 
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will  provide  an  animation  capability  for  viewing  aircraft  dynamic  motion  while  displaying  realtime  time-histoiy  data  of 
selected  dynamic  parameters.  The  resulting  coefficient  database  will  be  suitable  for  linear  analysis,  simulation,  or  any  of 
a  number  of  post-processing  analyses. 
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-  Rotational  Rates  ( P,Q,R ) 

-  Translational  Accelerations  (  a  and  P  ) 

-  Simulated  Tunnel  Test  Motion 
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-  Database  Entiy/Interpolation  with  (oCtrans  j  Ptrans)  variables 
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Each  module  will  provide  tabular  or  graphical  data  output.  This  provides  the  User  with 
insight  into  each  incremental  process  during  the  generation/ analysis  of  dynamic  data. 


SBIRN96-019 


Improved  Wnd  Tuimel  Data  Reduction  Procedure 


Progress  Report  #  3 


Sight,  Sound,  and  Motion 


SBIRN96-019 

(proved  Tiinnel  Data  Reduction  Procedure 


O  o 
o  d 


a 

o 

as 
B  si) 

OS 

Q  Z 


tTi 

rs 

S 

<N  ^ 

0  Vj 

X  4i 

SO  y 

it  I 

Os  Nv  ^ 
o  O)  ^ 

S  'C  ^ 

Ph  -R 

6-^  'P. 

QJ  <U 

s 

SI'S 

Sii 


R  'N 

Q  <: 


I 

ll  ^ 

1  s  -5 


325  East  Hillcrest  Dr.  #  100 
Thousand  Oaks,  CA  91360-5828 


SBIRN96-019 

iproved  Tunnel  Data  Reduction  Procedure 


I 


a 

o 


•G 

CJ 


o 

(/) 

bf) 

C/5 

V!) 

(U 

u 

a 

Gh 

o 

•a 

<U 

C/5 

G 

j 

CJ 

u 

u 

(/i 

a> 

c/5 

a 

§ 

HH 

u 

H 

o 

H 

o 

§ 

C/) 

g 

(U 

3 

’5 

’B 

2 

ii 

> 

(U 

bC 

s 

0 

s 

2 

1 

CA 

a 

a 

p:^ 

o 

o 

o 

a 

J? 

C/) 

H 

‘■P 

(j 

•^3 

o 

P 

G 

"§H 

HH 

’d 

g 

PQ 

c/3 

P^H 

o 

C/l 

C/) 

□  □  □  □  □ 


Sight,  Sound,  and  Motion 


SBIRN96-019 

mproved  Wind  Tunnel  Data  Reduction  Procedure 


Sight,  Sound,  and  Motion 


I 

s 

Ph 

c 

o 

•^3 

H 

S  ^ 
^  a 
52  Q 

0^  U 

HH  g 

^  s 

C/5  3 

H 

I 

’T3 

u 

s 

I 


a 

o 

*^. 

& 

•c 

C-) 

V) 

(S 

c« 

•  ^ 

o 

(U 

H 

’TS 

O 

’S 

a> 


a 

o 

*^3 

u 

>§ 


a 

o 

o 

s 

I 

a> 

W) 

a 

PlH 

"B 

c 

o 


a 

Q> 

li 

o 


^  P. 
2  2 
4-»  £j 

S 


U 


CJ 

o 

c/> 

bJD 

M 

W) 

o 

6^ 

<L| 

a 

a 

O 

1 

nd 

cn 

OJ 

d 

•  S 

b£  ^ 
d  1—1 

1 

<u 


biD 


a 


Sight,  Sound,  and  Motion 


VO 

I 


I 

u 

u 

s 

§ 

•o 

I 


Vi 


e 

I 


(i> 

V) 

a 

c/;) 

H 


o 

’S 


QJ> 

M-i 

o 

a 

.2 

i 

•I 

-a 

S 

PQ 


a 

J  si 
^  s-e 

ll'i 

«  i  ^ 

■s  si 

CO  (U  1^ 

o  S 

"O  ^ 

•s-y  ■ 


u 

I 

o 


^  -S 

^  r? 


\ 

11 

,4^ 


<L> 

.sq 

B 

nd 

o 

S 

B 


(L> 

23 

,  B 
^  6 
a  ji: 


O  ^ 

£1*'^ 

C3 

[2  '  Ch 

^  a  ° 

^  . 


a 


Oh  > 

3 

72 

^  O 

■“  I 

•l■■§ 
2  ^ 
S  O 


o 

Sh 

B 

O) 

.  I 


S-( 

R 

d) 


CO 

<D 


"d  oj 


Ul 

§ 

w 


CS 


i 
0-1 
d  ^ 

.  d  CO 


Sight,  Sound,  and  Motion 


SBIRN96-019 

proved  Wind  Tiumel  Data  Reduction  Procedure 


Sight,  Sound,  and  Motion 


Shown  here  is  the  prototype  GUI  for  entry  of  tunnel  test  and  flow  condition  data.  It  provides  the  Force- Moment  Synthesis 
module  with  input  parameters  for  generation  of  non-linear  aero  data.  The  color  plot  shows  the  input  waveform  for  oscillatoiy 
tests,  overlaid  with  ejqjected  goodness-of-fit  for  the  Fourier  approximation  given  S/N  ratio,  Cycles,  and  Samples/ Cycle. 


